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ABSTRACT 


In  this  report,  an  externally  pressurized 
axisymmetrical  foil  bearing  was  analyzed  in  order 
to  obtain  the  qualitative  effect  of  side  leakage 
in  a  foil  bearing.  One  or  two  circumferential 
pressure  sources,  symmetrically  placed,  were  used. 

In  the  analysis,  bending  stiffness  of  the  foil  was 
included  and  was  found  to  be  not  negligible.  Curves 
were  obtained  which  show  the  effects  of  the  various 
physical  parameters  on  the  gap  profile. 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  i 

TABLE  OF  CONTENTS  ii 

LIST  OF  ILLUSTRATIONS  Ill 

NOMENCLATURE  iv 

I.  INTRODUCTION  1 

USES  OF  FOIL  BEARINGS  1 

PROBLEMS  OF  ANALYSIS  1 

RESULTS  OF  PREVIOUS  STUDIES  2 

PURPOSE  OF  PRESENT  STUDY  2 

II.  ANALYSIS  AND  SOLUTIONS  3 

A.  FORMULATION  OF  THE  PROBLEM  FOR  ONE  SOURCE  3 

B.  APPROXIMATE  SOLUTIONS  6 

1.  Solution  for  Very  Thin  Foils  6 

2.  Solution  for  W2  <<  1  7 

C.  EXTENSION  TO  TWO  SOURCES  8 

III.  RESULTS  AND  CONCLUSIONS  10 

REFERENCES  26 

APPENDIX  I  1-1 

APPENDIX  II  II-l 


ii 


_ AMPEX 

LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  Axisymmetrical  Foil  Bearing  12 

2  Gap  Profiles  for  Perfectly  Flexible  Tape  13 

3  Variation  of  the  Indicial  Roots  as  Functions  of  G/a  14 

4  Tension  vs.  Flow  Rate  for  Source  at  b  =  0,3  15 

5  Tension  vs.  Flow  Rate  for  Source  at  b  *  0,7  16 

6  Tension  vs.  Flow  Rate  for  Source  at  b  =0.85  17 

7  Tension  vs.  Flow  Rate  for  Source  at  b,  =  0.95  18 

8  Gap  Profiles  for  Source  atbj  =  0.3  19 

9  Gap  Profiles  for  Source  atbx  =  0.7  20 

10  Gap  Profiles  for  Source  at  bj  =  0.85  21 

% 

11  Gap  Profiles  for  Source  at  bj  =  0.95  22 

12  Typical  Gap  Profiles, bj  =  0,85  23 

13  Typical  Gap  Profiles ,bj  =0.3  24 

14  Typical  Gap  Profiles,  bj  =  0.7  25 


iii 


NOMENCLATURE 


b  half-length  of  the  foil  cylinder  (in.) 

b'  distance  of  source  from  foil  center  line  (in.) 

b  b '  /b 

c  unstretched  circumference  of  foil  (in.) 

h  gap  width  between  the  bearing  surfaces  (in.) 

hQ  gap  width  under  foil  center  line  (in.) 

p  absolute  pressure  in  lubricating  air  film  (psi) 

pa  ambient  pressure  (psi) 

p  absolute  pressure  at  source  (psi) 

r  radius  of  foil  cylinder  (in.) 

t  foil  thickness  (in.) 

w  deflection  of  foil  from  gap  hD  (in.) 

x  axial  coordinate  measured  from  foil  center  line  (in.) 


A 

B 

C 


flexibility  parameter. 


12  ( 1-v2  )b‘t 
R2t2 


edge  curvature  parameter,  t/hD 


dimensionless  parameter. 


PaR2 

Et 


C  -C4  constants  of  integration 

Et  ^ 

D  bending  modulus,  - - -  (lb-in.) 

12 ( 1- v  2 ) 

E  modulus  of  elasticity  (psi) 


F 


dimensionless  tension  parameter, 


V L' 

h0Et 


iv 


a  /A/ 4' 

6  ctbj 

n  x/a 

n j  — n 5  possible  real  and  imaginary  parts  of  ri 

X  index  needed  for  solving  differential  Equation  (26) 

X^  possible  values  for  \ 

u  viscosity  of  air  ( lb-sec/in, 2 ) 

v  Poisson’s  ratio  for  foil 

P„  density  of  air  at  ambient  pressure  (lb  mass/in,3) 

8. 
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I.  INTRODUCTION 


A  bearing  with  at  least  one  wall  made  of  a  thin  flexible  material 
is  known  as  a  foil  bearing.  Actually,  in  all  bearings,  both  surfaces 
are  flexible  to  some  extent.  Thick-walled  bearings  have  been  analyzed 
by  assuming  rigid  surfaces;  and  thin-walled  bearings  have  been  analyzed 
by  assuming  perfectly  flexible  surfaces  (i.e.,  surfaces  having  no 
bending  stiffness).  These  two  limiting  cases  are  equivalent  to  the 
limiting  cases  of  infinite  thickness  and  zero  thickness  of  the  bearing 
wall.  In  this  report,  the  thickness  will  be  considered  as  small  but 
not  negligible. 

USES  OP  FOIL  BEARINGS 

In  the  most  general  sense,  a  foil  bearing  is  used  whenever  a  moving 
foil  or  thin  sheet  (commonly  called  a  "web"  in  the  paper  industry)  has 
to  be  supported.  Because  of  the  light  weight  of  most  foils,  a  simple 
stationary  support  is  often  adequate.  If,  in  addition  to  support,  a 
change  in  the  direction  of  motion  of  the  foil  is  desired,  the  friction 
of  a  stationary  support  is  likely  to  be  too  high.  In  this  case,  a 
support  that  moves  with  the  foil  (rotating  support)  or  a  lubricating 
film  can  be  used.  The  stationary  and  rotating  types  of  supports, 
without  lubrication  between  the  foil  and  the  support,  have  been  widely 
used  in  the  paper,  textile,  plastics,  and  metal  industries,  and  also 
on  magnetic  information  storage  devices  such  as  tape  recorders.  To  the 
author's  knowledge,  lubricating  films  between  foil  and  stationary  guides 
have  only  been  used  in  tape  recorders,  metal  processing!]!] ,  and  in  test 
equipment  for  loading  rotating  shafts[2],  In  all  these  cases,  the 
equipment  operates  in  air,  and  atmospheric  air  is  therefore  the  lub¬ 
ricant  . 

Foil  bearings  can  also  be  used  to  support  entire  rotors,  replacing 
either  ball  or  rigid  surface  fluid  bearings  for  this  purpose.  The  ad¬ 
vantages  of  foil  bearings  are  that  they  use  fluid  films  for  developing 
forces  and  therefore  have  the  low  friction,  low  noise,  and  long  life 
of  rigid  surface  fluid  film  bearings;  but,  like  ball  bearings,  they  do 
not  have  the  self-excited,  whirl  problem  common  to  rigid  surface  fluid 
film  bearings.  It  is  certain  that  as  the  behavior  of  foil  bearings 
becomes  better  understood  and  design  tools  other  than  trial  and  error 
become  available,  many  additional  uses  for  these  bearings  will  be  found. 

PROBLEMS  OF  ANALYSIS 

Analysis  of  foil  bearings  involves  simultaneous  solution  of  the 
Reynolds  equation  (relating  the  pressure  in  the  fluid  film  to  the 
parameters  of  that  film)  with  the  foil  equilibrium  equations.  The  re¬ 
sulting  equation  is  a  sixth-order,  nonlinear,  nonhomogeneous  differential 
equation,  for  which  no  solution  has  been  found.  Because  of  the  high 
order  of  the  equation,  finite  difference  solutions  are  difficult,  and 
emphasis  has  therefore  been  on  approximate  solutions.  However,  by 
treating  the  problem  as  a  one-dimensional  problem  and  assuming  the  foil 
to  be  a  perfectly  flexible  membrane,  the  differential  equation  can  be 
reduced  to  a  third-order  nonlinear  differential  equation. 
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It  should  be  noted  that  even  the  rigid  surface  problem  is  diffi¬ 
cult  when  the  lubricating  film  is  a  gas.  The  state-of-the-art  of  foil 
bearings  is  approximately  that  of  gas  bearings  seven  years  ago;  i.e., 
some  approximate  solutions  have  been  obtained,  and  some  understanding 
of  their  behavior  has  been  gained. 

RESULTS  OP  PREVIOUS  STUDIES 

Because  of  the  difficulty  of  analysis  and  the  limited  number  of 
applications,  few  studies  have  been  made  of  foil  bearings,  and  all  of 
these  have  been  of  self-acting  bearings.  Blok  and  Van  RossumC3l, 
who  did  a  perturbation  study  of  oil-lubricated  foil  bearings,  found 
that  the  gap  and  the  pressure  were  approximately  constant.  However, 
their  experiments  showed  that  cavitation  occurred  in  the  exit  flow 
region.  In  their  study,  as  well  as  in  perturbation  analyses  of  foil 
bearings  by  Langlois  and  Kohl!  ^3  and  Grossly],  the  foil  was  assumed  to 
be  perfectly  flexible  and  infinitely  long  in  the  axial  direction. 

These  two  assumptions  greatly  simplified  the  foil  equilibrium  equation, 
and  also  permitted  axial  flow  and  axial  variations  in  gap  and  pressure 
to  be  neglected.  However,  the  foil  cannot  be  considered  perfectly 
flexible  whenever  the  second  derivative  of  curvature  is  of  the  same 
order  as  the  ratio  of  pressure  to  bending  stiffness.  This  condition 
will  occur  in  the  entrance  and  exit  flow  regions, 

PURPOSE  OP  THE  PRESENT  STUDY 

In  a  finite  bearing,  some  of  the  flow  is  in  the  axial  direction 
and  is  lost  at  the  side  edges.  This  side  leakage  causes  the  pressure 
and  the  gap  to  decrease  from  that  predicted  for  infinitely  wide  bearings. 
In  the  present  study,  this  problem  was  formulated  in  order  to  obtain  a 
qualitative  idea  of  the  effect  of  side  flow  and  finite  foil  stiffness. 

The  results  of  this  study  can  also  be  applied  to  externally  pressurized 
foil  bearings  in  regions  with  primarily  axial  flow  and  almost  constant 
circumferential  curvature. 

In  order  to  eliminate  flow  and  variations  in  gap  and  pressure  in 
the  circumferential  direction,  the  foil  is  assumed  to  be  wrapped  com¬ 
pletely  around  a  cylinder,  with  a  constant  pressure  circumferential 
line  source  placed  under  the  foil  center  line  in  order  to  provide  axial 
flow  (Figure  1,  with  b'=  0),  The  pressure  drop  to  ambient  at  the  edge 
is  expected  to  cause  a  similar  decrease  in  gap  from  the  center  to  the 
edge.  Near  the  edges,  the  bending  stiffness  of  the  foil  will  cause  it 
to  be  concave  upwards.  This  anticlastic  curvature  effect  is  always 
found  in  the  bending  of  thin  plates!!?]. 
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II.  ANALYSIS  AND  SOLUTIONS 


A.  FORMATION  OF  THE  PROBLEM  FOR  ONE  SOURCE 


A  line  source  of  strength  Pj  atmospheres  is  placed  around  a  cyl¬ 
inder  of  radius  R.  A  thin  foil  of  initial  length  c,  width  2b,  and 
thickness  t,  is  wrapped  around  the  cylinder  so  that  the  line  source 
is  under  the  center  line  of  the  foil  (Figure  1,  with  b*  =0). 


The  gap  h0  under  the  center  line  is  related  to  the  circumferential 
tension  T0  along  the  center  line  by  the  formula. 


c 

2  if 


1  +  _o 
Et 


-  R, 


(1) 


where  E  is  Young's  modulus  for  the  foil.  The  foil  is  assumed  to  be 
made  of  an  isotropic  material  with  Poisson's  ratio  v.  Since  the  foil 
is  thin  and  the  deflections  from  the  circular  cylindrical  shape  are  of 
the  same  order  as  the  thickness,  we  can  use  the  plane  stress  approxi¬ 
mations  to  the  elastic  equations.  Under  these  assumptions,  the  equi¬ 
librium  equations  for  the  foil  are  the  same  as  those  for  symmetrical 
cylindrical  shellsC8],  Thus, 


Et3  ,  d\/_  +Et 
12 ( 1-v2 )  *  dx4  R2 


^  -  (P-Pa) 


(2) 


dNx 

dx 


qx  = 


dw 
P  dx 


(2a) 


where  w  is  the  deviation  from  the  cylinder  r  =  R  +  h0,  and  where  Nx 
is  the  tangential  stress  resultant  and  qx  is  the  component  of  dis¬ 
tributed  external  loading  in  the  x-direction.  The  loading  qx  is  small 
except  near  the  edge,  where  the  slope  can  be  large,  but  Nx  must  decay 
to  zero  at  the  edges,  which  are  unconstrained.  Thus,  the  normal  com¬ 
ponent  of  Nx  can  be  neglected,  as  was  done  in  Equation  (2).  Equation 
(2a)  can  be  used  to  find  Nx,  once  qx  has  been  determined. 

Since  the  problem  has  axial  symmetry,  we  are  concerned  with  one¬ 
dimensional  flow.  The  governing  equation  of  the  air  film  will  be  the 
Reynolds  equation  for  two  surfaces  that  are  stationary  with  respect 
to  each  other.  Thus, 


d_ 

dx 


d£ 

dx. 


0 


(3) 
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Equations  (2)  and  (3)>  together  with  appropriate  boundary  conditions, 
are  sufficient  for  finding  the  gap  and  pressure  profiles. 

The  boundary  conditions  due  to  symmetry  about  the  foil  center  line 
are ,  at  x  =  0 : 


dw  _ 
dx 


0, 


V  =  ^5  =  0 
dx3 


(4) 


The  pressure  and  free  edge  conditions  give  us,  at  x  =  b: 


d2w  _ 
dx2 


0, 


and 


Et3 

12 ( 1-v2 ) 


*  +  Ew  =  T. 

dx4  R2  R 


At  x  =  0, 


Et3 

12 ( 1-v2 ) 


d4w  Et 
dx4  R2 


_a  - 

R 


Pa(pi"  D 


(5) 


(6) 


(7) 


The  conditions  of  Equation  (5)  result  from  the  vanishing  of  the  moment 
and  shear  resultants  at  the  edge,  and  those  of  Equations  (6)  and  (7) 
result  from  the  given  pressures,  Pj  (p,/pa)  atmospheres  at  the  center 
line,  and  ambient  pressure  pa  at  the  edge. 

Making  the  substitutions 


H 


X  =  ~,  and  P  =  — 
b  Pa 


(7a) 


in  Equations  (2)  through  (7),  we  obtain 


(h3p-p'  )'  =  0 


(8) 


WiV(0)  +  AW( 0 )  =  A  [F-C(P  -1)] 


(13) 


In  the  above  equations,  the  prime  represents  differentiation  with 
respect  to  X,  and  A,  B,  F,  and  C  are  given  by 


a  .  HI1!?2!*4.,  B  - 

R2t2  hn 


T_R 
F  =  0 


p  R2 

C  = 


In  the  present  form,  Equations  (8)  and  (9)  are  difficult  to  solve 
unless  we  perturb  about  the  incompressible  solution.  If  Pj  Is  of 
order  1  +  6,  where  6  <<  1,  then  a  one-term  perturbation  is  satisfactory. 
Since  in  most  self-acting  foil  bearings,  T/paR  is  much  less  than  one, 
the  perturbation  method  can  be  used,  keeping  only  the  lowest  order 
terms,  which  correspond  to  the  incompressible  solution. 

If  the  lubricating  film  Is  incompressible,  Equation  (8)  reduces  to 


(15) 


where  -  G/C  is  the  constant  coming  from  the  first  integration  of 
Equation  (8)  and  is  related  to  the  flow  rate  Q  by 


G  =  12ub 
^  Pa^o3Pa 


(16) 
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Equations  (9)  and  (15)  can  be  combined,  yielding 


Wv  +  AW/ 


GA 

(l-W)3 


(17) 


B.  APPROXIMATE  SOLUTIONS 

Two  approximate  solutions  are  of  interest:  (1)  for  very  thin  foils, 
and  (2)  for  W2  <<  1. 

1.  Solution  for  Very  Thin  Foils 

When  t  is  small,  A  is  large,  and  Equation  (17)  can  be  reduced 

to : 


W 


/ 


G 

"(1-W)S 


(18) 


Integration  of  Equation  (18)  gives 


(1-W)4  -  (1-We)4  =  4G  (1-X), 


(19) 


where  We  is  the  deflection  at  X  =  1.  The  boundary  conditions  are:. 

W(0)  =  -  C(P:-1)  +  F  (20) 

and 

W(l)  =  We  =  F  (21) 


By  definition  of  W  and  hQ,  we  have 
W(0)  =  0 


(22) 


Evaluating  Equation  (19)  at  X  =  0  gives 


1  _  ( 1-We ) 4  =  4g 
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(23) 


This  solution  is  expected  to  be  valid  everywhere  except  at  the 
center  line,  X  =  0,  and  near  the  edge,  X  =  1.  In  those  two  regions 
W1  =  0(AW),  and  the  conditions  for  the  derivatives  of  W  are  not 
satisfied.  Gap  profiles  for  several  values  of  We  are  shown  in  Figure  2. 

Equations  (18)  and  (22)  give,  at  X  =  0+, 


/ 

W 


=  G 


(24) 


Since  the  flow  is  outward  from  the  source,  the  gap  is  symmetrical 
about  the  W-axis.  Thus  for  X  =  0”, 


W  '  =  -  G 


(25) 


That  is,  the  perfectly  flexible  foil  has  a  crease  along  its  center 
line.  Since  the  perfectly  flexible  foil  has  no  bending  stiffness,  the 
anticlastic  curvature  does  not  occur  at  the  edge  (compare  with  Figures 
8  through  14), 

* 

2.  Solution  for  W2<<  1 

When  W2<<  1,  the  right  hand  side  of  Equation  (17)  can  be  expanded, 
giving 


WV  +  AW  GA  (1+3W) 


(26) 


The  solution  of  Equation  (18)  has  the  form 


W 


(  5 

V 

L- 

K=1 


CkexkaX 

j 


1 

3 


(27) 


where  the  X^' s  represent  the  roots  of 


Xs  +  4X  -  ill  =  0, 

a 


(28) 


and 


a  = 


(29) 
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Since  Equation  (28)  has  one  real  root  (see  Figure  3),  Equation  (27) 
can  be  rewritten  as  follows: 


W  =  — i  +  K1er,1°<X  +  en2aX(K2sinn3aX+K3cosn3aX) 


+  en4“X(K4sinnsaX+K5cosri5aX) 


(30) 


where  the  ni's  are  real  and  have  been  formed  from  the  X^'s. 

Figure  3  shows  the  n^'s  as  functions  of  G/«  .  The  boundary 
conditions,  Equations  (10)  through  (13)  and  Equation  (22),  determine 
Kj  through  K5  and  G,  respectively.  However,  we  have  seven  equations 
and  six  unknowns,  allowing  us  to  use  one  equation,  Equation  (13),  as 
a  relation  between  T0  and  P,  .  Since  the  n^'s  are  dependent  upon  G, 
the  boundary  equations  are  transcendental.  For  this  reason  in  the  com¬ 
puter  programs  for  this  problem,  G  was  assumed  to  be  known,  and  T0  and 
P,  were  calculated.  (For  calculations  see  Appendix  I.)  Gap  pro¬ 
files  have  been  plotted  only  for  the  two-source  problem,  of  which  this 
problem  can  be  considered  a  special  case. 

C.  EXTENSION  TO  TWO  SOURCES 

As  an  extension  of  the  above  problem,  two  circumferential  line 
sources,  located  at  x  =  ±b',  were  used.  There  is  no  flow  between 
these  sources,  and  the  equation 


W1V  +  AW  =  A[F-C(Pj -1) ]  (31) 

holds  in  the  edge  regions.  For  x  > bj  ,  we  have  X  >  b, ,  where  bj  =  b'/b t 
and  Equation  (17)  is  the  governing  differential  equation  for  W.  In 
addition  to  the  boundary  conditions  of  Equations  (10),  (11)  and  (22), 
we  have  the  following  continuity  conditions  at  X  =  b3: 

W(b1”)  =W(b1+),  W'(b:“)  =U'(b1+)t 

w'^bj-)  =w"(b1  +  ),  w'"  (bj")  -  Ww  (bj  +  )  (32) 

Equation  (12)  is  replaced  by 

WiV(b1+)  +  AW(bJ+)  =  A[F-C(P  -1)]  (33) 


8 


The  approximate  solution  developed  in  Equations  (26)  through  (30) 
can  be  used  in  the  region  X  >  bj  .  The  equations  for  the  constants  of 
integration  are  developed  in  Appendix  I.  The  curves  relating  T  to  G, 
with  parameters  B  and  A,  are  shown  in  Figures  4  through  7,  each  figure 
representing  a  different  source  position.  Some  typical  profiles  are 
shown  in  Figures  8  through  14, 

For  more  than  two  sources,  a  numerical  integration  procedure  is 
necessary,  A  computer  program  for  numerical  integration  of  Equation 
(17)  for  an  arbitrary  number  of  sources  has  been  developed.  By  using 
this  program  we  can  adjust  the  number  of  sources  and  their  relative 
strengths  to  give  any  desired  flow  rate  profile. 
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III.  RESULTS  AND  CONCLUSIONS 


If  the  foil  were  perfectly  flexible  (A-*-®),  the  foil  tension  T0 
would  be  related  to  P2  by 


p  (^-1)  =  -i 
a  R 


The  bending  stiffness  of  the  foil  provides  some  normal  shear,  giving 
us  Equation  (7).  However,  the  left  hand  side  of  Equation  (7)  will  be 
small,  so  that  P2  and  T0  are  related  by  Equation  (34),  except  for 
higher  order  terms.  Thus  of  P2,  T0,  h0,  and  Q,  we  can  specify  h„  and 
one  of  the  others,  and  can  solve  for  the  remaining  variables.  Given 
Pj  and  h0,  or  T0  and  hOJ  we  can  also  adjust  Q  by  using  two  sources  and 
varying  their  distance  from  the  foil  edges.  Within  the  limitations 
mentioned,  Pj  and  T0  can  be  adjusted  to  determine  h0.  As  the  tape 
stiffness  increases,  the  edge  curling  effect  becomes  more  pronounced, 
as  is  shown  in  Figures  12,  13,  and  14.  In  Figures  12  and  13,  the 
effect  of  increasing  thickness  is  shown  for  two  source  positions; 
in  Figure  14,  the  modulus  of  elasticity  is  varied. 

The  flexibility  of  the  foil  is  determined  by  the  relative  impor¬ 
tance  of  the  term 


Et J 

12 ( 1-v  2 ) 


compared  to  the  terms 


(36) 


in  Equation  (2).  As  the  term  (35)  becomes  negligible,  the  foil  be¬ 
comes  a  membrane.  Equation  (2)  can  be  transformed  into  the  dimension¬ 
less  form  given  by  Equation  (9).  The  flexibility  parameter  in  Equation 
(9)  is  represented  by  A,  As  A  increases,  the  foil  becomes  more  flexible. 
Another  Influence  on  the  flexibility  is  the  edge  curvature  parameter  B. 
Decreasing  B  increases  the  foil  flexibility  by  decreasing  the  curvature 
change.  The  bending  term  becomes  important  if  A.  becomes  small,  if  the 
higher  derivatives  of  W  are  large  compared  to  W,  or  if  a  combination 
of  these  two  effects  occurs.  Examining  A  (Equation  14),  we  see  that 
decreasing  the  thickness,  decreasing  the  radius  of  the  cylinder,  or 
increasing  the  foil  width  will  make  the  foil  more  flexible.  Decreasing 
the  thickness  or  increasing  the  central  gap  also  decreases  B  (Equation 
14),  again  making  the  foil  more  flexible.  A  good  comparison  can  be 
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AMPEX 


made  using  Figures  4  through  11.  First,  F,  B,  and  A  must  be  calcu¬ 
lated  from  the  physical  parameters  - (F  and  C[P1-1]  are  sufficiently 
close  to  allow  one  curve  to  represent  both  parameters.)  Using  each 
set  of  curves  corresponding  to  a  given  source  location  (Figures  4 
through  7),  we  find  the  proper  value  for  G.  The  gap  profile  can  then 
be  found  using  Figures  8  through  11,  with  interpolation  for  G.  Or,  it 
can  be  calculated  using  Equations  (la)  and  (30),  with  Cj,  C2,  and 

through  K5  given  by  Equations  (3a),  (4a),  and  (6a)  through  (10a)  in 
Appendix  I,  where  the  n^'s  can  be  calculated  by  the  method  shown  in 
Appendix  II. 

As  the  sources  are  moved  nearer  the  edges,  we  find  from  Figures  8 
through  11  that  the  effect  of  flow  rate  on  the  gap  profile  diminishes. 
This  result  is  due  to  the  increased  effect  of  foil  stiffness  near  the 
edges,  since  the  gap  tends  to  remain  nearly  constant  between  the 
sources  and  resists  rapid  changes  outside  that  region.  The  central 
gap  hQ  can  be  varied  by  varying  TQ  and  P1§  Increasing  T0  tends  to 
decrease  hQ,  while  increasing  Pj  tends  to  increase  h0. 

Since,  in  a  self-acting  foil  bearing,  the  side  flow  varies  from 
zero  under  the  foil  center  line  to  a  maximum  at  the  edges,  using  more 
than  one  source  allows  us  to  improve  our  approximation  to  the  side 
flow.  The  major  limitation  on  the  approximate  solution  presented  is 
that  for  the  assumption  W2  <<  1  to  be  valid,  the  minimum-to-central 
gap  ratio  cannot  be  greater  than  0.5.  The  solutions  presented  are 
useful  in  any  externally  pressurized  foil  bearing  which  has  an  almost 
constant  radius  and  little  circumferential  flow. 
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±i>?3.,ri4  ±  iT»5 
ARE  ROOTS  OF  r?°  +  4  n  -12  G/«C=0. 


Figure  3..  VARIATION  OF  THE  INDICIAL  ROOTS  AS  FUNCTIONS  OF  G/*C 
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Figure  6.  TENSION  VS.  FLOWRATE  FOR  SOURCE  AT  b,  =0.85 
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Figure  8.  GAP  PROFILE  FOR  SOURCE  AT  b  =  0.3 
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Figure  12.  TYPICAL  GAP  PROFILES,  b  =  0.  85 
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Figure  13,  TYPICAL  GAP  PROFILES, b 
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APPENDIX  I 


Equation  (29)  can  be  integrated  to  give 

W  =  F  -  CCPj-l)  +  Cjcosh  aX  cos  aX  +  C2sinh  ax  sin  “X 

+  C3sinh  aX  cos  aX  +  C^cosh  aX  sin  aX  (la) 

Substituting  into  Equation  (10),  we  find  that 

C3  =  C4  =  0  (2a) 

Equation  (la)  holds  for  X  i  bj .  Outside  this  region.  Equation  (30) 
will  be  used  as  an  approximate  solution  of  Equation  (17).  Substi¬ 
tuting  Equation  (17)  into  (11),  letting  abj  =  3,  gives 

Kj  nj  2  en  1  a+en2a  {  K2  [  ( n 2 2  — n  3  2  )  sin  n3a+2n2n,  cos  n3<*] 

+K3  C (n22— n32  >  cos  n3a-  n2 n 3  sin  n3«]  } 

+eri4a  {k4  [(n4z-n52)  sin  n5a  +2n4n5  cos  n5a] 

+K5  [(n42-n52)  cos  n5a  -2n4n5  sin  n5a]  } 


Ki  n  x  3 en  1  “+^2°  (k2  [n2  (n22-3n32)  sin  n3a 

+n3  (3n22_T132  )  cos  ^ 3 a  +k3  C^2  (n22-3n32  )  cos  n3« 

-n3  (3n22-n32  )  sin  ^a]}  +eri401  {  K4  [n4(r>42-3n5  2)  sin  n,.a 

+n5(3\2-n53)  cos  n5a  +Kg [n4  ( \2-3n52  )  cos  n5a 

-n5 (3n42-n52  )  sin  n5ct]}  =  0  ( ) 

Substituting  Equations  (la)  and  (28)  into  Equations  (22)  and  (30) 
gives 

c  +  P  -  C  (Pj-1)  =  o  (5a) 

-  — +K. en 1 S+en2 ^ (K  sin  n  B+K  cos  n  B) 

^  I  2  3  3  3 

+en|*^(K  sin  n  B+K,  cos  n  3) 

4  5  5  5 

=  c  (cosh  6  cos  8-1)  +  C2  sinhB  sin  B,  (6a) 


1-2 


eT1ie+en23[K2(n2  sin  n3B+n3  cos  n3 6) 

+K  (n  cos  n  8-n  sin  n  B)] 

3  2.  J  3  ^ 

+eni*6[K4(n4  sin  n,.B+n5  cos  nsB)  +Ks(n4  cos  n5B-n5  sin  ngS)] 
=  C, ( sinh  8  cos  8  -cosh  8  sin  8)  +C, (cosh  B+sinh  8  cos  e), 

X  *- 

(7a) 


Kifij26"1!  6+en2e  {K2[(n22-n32)  sin  n3B+2n2n3  cos  n  3  B  3 
+K3[(n22-n32)  cos  n3B-2n2n3  sin  n 3 B ] } 

+en4e{K4  [(n42-n52  )  sin  n5B+2n4ri5  cos  n5B] 

+K5[(n42-ri52)  cos  n5B-2n4n5  sin  n 5 B 3 } 

=  2(Cz  cosh  B  cos  B-Cj  sinh  6  sin  B),  (8a) 
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Kiri13eniS+en2S  {K2[n2(n22-3n32)  sin  n3  e+n3  ( 3n22 -n32  ) 
cos  n 3 e  3  +K3 [n2 (n22-3n32 )  cos  n 3 e -n 3 (3r>22-r'32 )  Sin  n 3 B ] } 
+en*+3  {K4[n4(nl+2-3n52)  sin  n5B+n5(3n42-n52 )  cos  n 5 B 3 


+K 5 [ rm  (  ni+2-3n  52  )  cos  n5B-n5(3n42-n52  )  sin  n 5 3 3  } 

=  2  [-C  (sinh  B  cos  B  +  cosh  B  sin  B) 

+C2(sinh  B  cos  B  -  cosh  B  sin  B)]  (9a) 

K1ni‘leni8+en2S  {K2[(n2‘*-6n22n32+n34)  sin  n3B 
+  1*  n2n3(  n22-n32)  cos  n3B]  +K  3C  (  n24-6  n22n 32+ n 34)  COs  n3B 
n2  n3  (  n22- n32  )  sin  n3B]}  +eni<6  {k4 [ ( n44-6n42 n52+n 54 ) 
sin  n5B  +  4n4n5(n42-n52)  cos  n 5 3 ]  +K5[ ( n44_6n42n52+n54) 
cos  n5B-4n4n5(n42-n52)  sin  n 5 3] } 

=  — 4 ( C  j  cosh  B  cos  B  +  C2  sinh  B  sin  B)  (10a) 

■Equations  (3a),  (4a),  and  (6a)  through  (10a)  form  a  set  of  seven 
simultaneous  linear  algebraic  equations  for  Cj,  C2,  and  Kj  through  K5. 
Equation  (5a)  gives  the  relation  between  F  and  Plt  The  above  equations 
were  solved  for  various  values  of  A,  B,  and  G,  using  the  7090  computer 
at  David  Taylor  Model  Basin. 
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APPENDIX  II 


The  indicial  equation, 


X5  +  A  -  3AG  =  0 


(lb) 


cannot  be  solved  by  algebraic  means.  Rolle's  theorem  gives  us  a 
restriction  on  the  location  of  the  roots.  Defining  K  by 


K  =  X5  +  AX  -  3AG 


(2b) 


the  maxima  and  minima  of  K  are  given  by 


—  =  5X4  +  A  =  0  (3b) 

dx 

Since  Equation  (3b)  has  no  real  roots,  K  is  a  monotonically  in¬ 
creasing  function  and  must  vanish  at  one  point.  For  x  =  0,  K  is 
negative.  Therefore,  the  real  root  of  Equation  (lb)  is  positive. 
Calling  this  root  Xj  =  n x a ,  and  setting  X  =  na,  we  can  write 
Equation  (lb)  as 

(ti-n1)(n4+C1n3+C2n2+C3n+C4)  =  n5+4n-12  £  =  0  (4b) 


where 


Ci  =  ni>  C2  =  ^l2’  C3  =  ni3»  and  C4  =  ni4+il  (5b) 


n 

If  we  set  u  =  “  ,  the  equation  for  the  complex  roots  becomes 
l 

u4  +  u3  +  u2  +  u  +  1  +  -i—  =  0  (6b) 

V 

The  value  for  rij  can  be  found  from  Equation  (2b)  by  graphical 
means  and  can  be  refined  by  using  the  formula 


II-l 


(7b) 


AM: 


4nT  5+12— 
1  a 

5n74+4 


where^  is  the  graphical  solution.  Equation  (6b)  can  be  solved 
exactly[6] . 
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Attn:  Librarian  1 

Mr.  William  D.  Stimmel 
Research  Laboratories  Library 
Allis-Chalmers  Mfg.  Co. 

Milwaukee  1,  Wisconsin  1 

J.  W.  Pecker  Division 
American  Optical  Company 
4709  Baum  Blvd. 

Pittsburgh  13,  Pennsylvania  1 

American  Society  of 

Lubrication  Engineers 
5  North  Wabash  Avenue 
Chicago  2,  Illinois  1 

Chairman 

Research  Committee  on  Lubrication 
The  American  Society  of 
Mechanical  Engineers 
United  Engineering  Center 
345  East  47th  Street 
New  York  17,  New  York  2 

Mr.  James  R.  Kerr,  President 

Lycoming  Division 

AVCO 

Stratford,  Connecticut  1 

Research  Precision  Mechanisms 
Division  of  the  Barden  Corp, 

4  Old  Newtown  Road 
Danbury,  Connecticut 
Attn:  B.  L.  Mims 

Vice  President-Gen,  Mgr.  1 

Beemer  Engineering  Company 
Industrial  Park 

Fort  Washington,  Pennsylvania  1 

Utica  Division 

The  Bendix  Corporation 

211  Seward  Avenue 

Utica,  New  York 

Attn:  Mr.  Russell  T.  DeMuth 

Supervisory  Engineer  1 


Bendix  Aviation  Corporation 
Research  Laboratories  Division 
Bouthfield,  Michigan 
Attn:  Mr.  Ralph  H.  Larson  1 

Mr.  C.  R.  Adams 
Physics  Technology  Department 
Aero-Space  Division 
The  Boeing  Company 

Seattle  24,  Washington  1 

Bryant  Chucking  Grinder  Company 
60  Clinton  Avenue 
Springfield,  Vermont 

Attn:  Mr.  Raold  Cann  1 

Cadillac  Gage  Company 
P.O.  Box  3806 
Detroit  5,  Michigan 

Attn:  Mr.  J.  Taylor,  Proj .  Mgr.  1 

Mr.  Dewey  J.  Sandell 
Director  of  Development 
Carrier  Res.  &  Dev.  Company 
Carrier  Parkway 

Syracuse,  New  York  1 

Chance  Vought  Corporation 
P.O.  Box  5907 
Dallas,  Texas 
Attn:  Mr.  R.C.  Blaylock 

Vice  Pres . (Engineering)  1 

Chrysler  Corporation 
Defense  Operations 
P.O.  Box  757 
Detroit  31,  Michigan 

Attn:  Mr.  C.  W.  Snider  1 

J,  D.  Mamarchev  &  Associates 
Consulting  Engineers 
3908  Main  Street 

Houston  2,  Texas  1 

Mr.  B.  W .  Birmingham 
Cryogenic  Engineering  Laboratory 
National  Bureau  of  Standards 
Boulder,  Colorado  2 

Curtiss  Wright  Corporation 
Wright  Aeronautical  Division 
Department  8332 
Wood  Ridge,  New  Jersey 
Attn:  W.  J,  Derner 

Chief  Project  Engineer  1 
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Daystrom  Pacific 
9320  Lincoln  Boulevard 
Los  Angeles  45,  California 
Attn:  Robert  H.  Smith 

Special  Projects  Engineer  1 

Ford  Instrument  Company 

31-10  Thomson  Avenue 

Long  Island  City  l,,New  York 

Attn:  Mr.  Jarvis  1 

Mr.  Adolf  Egli 
Ford  Motor  Company 
Engineering  &  Research  Staff 
P.O.  Box  2053 

Dearborn,  Michigan  1 

Dr.  John  E.  Mayer,  Jr. 

Non-Metallic  Section 
Applied  Science  Department 
Scientific  Laboratory 
Ford  Motor  Company 
P.O.  Box  2053 

Dearborn,  Michigan  1 

AiResearch  Manufacturing  Div. 

The  Garrett  Corporation 
9851  S.  Sepulveda  Boulevard 
Los  Angeles,  California 
Attn:  Jerry  Glaser,  Supervisor 

Mechanical  Lab,  Dept. 93-17  1 

General  Atomics  Division 
General  Dynamics  Corporation 
P.O.  Box  608 

San  Diego  12,  California 

Attn:  Mr.  F.  W.  Simpson  1 

Bearing  and  Lubricant  Center 

General  Engineering  Laboratory 

General  Electric  Company 

1  River  Road 

Schenectady,  New  York 

Attn:  G.  R.  Fox,  Manager  2 

Mr.  L.  W.  Winn 

General  Electric  Company 

Aircraft  Accessory  Turbine  Dept. 

950  Western  Avenue,  Bldg.  3-74 
Lynn,  Massachusetts  1 


Research  Laboratories 
General  Motors  Corporation 
General  Motors  Technical  Center 
12  Mile  and  Mound  Roads 
Warren,  Michigan 
Attn:  Mr.  E.  Roland  Maki 

Mechanical  Dev.  Dept,  1 

A.  C.  Spark  Plug  Division 

General  Motors  Corporation 

Milwaukee  1,  Wisconsin 

Attn:  Allen  Knudsen  1 

Mr.  Walter  Carow 
Kearfott  Division 
General  Precision  Inc. 

1150  McBride  Avenue 

Little  Falls,  New  Jersey  1 

Gruman  Aircraft  Engineering  Corp. 
Bethpage,  Long  Island,  New  York 
Attn:  Mr.  David  W.  Draig,  Jr. 

Mechanical  Design  Section 
Engineering  Department  1 

Hydronautics ,  Incorporated 
Pindell  School  Road 
Howard  County 

Laurel,  Maryland  1 

International  Business  Machines 

Corporation,  Research  Laboratory 
San  Jose,  California 

Attn:  Dr.  W.  E.  Langlois  2 

Mr.  B.  A.  Napier,  Engr.  Div.  Mgr. 
Lear,  Incorporated 
110  Ionia  Avenue  N.W. 

Grand  Radids,  Michigan  1 

Mr.  L.  R.  Barr,  Vice  Pres. 

Lear-Romec  Division 
Abbe  Road 

Elyria,  Ohio  1 

Dr.  Calus  G.  Goetzel,  D/53-30 
Bldg.  201,  Plant  2,  Palo  Alto 
Lockheed  Missiles  &  Space  Co. 

P.O.  Box  504 

Sunnyvale,  California  1 
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Dr.  J.  S.  Ausman 
Litton  Systems,  Inc. 

5500  Canoga  Avenue 

Woodland  Hills,  California  2 

Mr.  Don  Moors 
Litton  Systems,  Inc. 

5500  Canoga  Avenue 

Woodland  Hills,  California  1 

Mr.  A.  N,  Thomas 

Astro  Division 

The  Marquardt  Corporation 

16555  Saticoy  Street 

Van  Nuys,  California  1 

Mr,  Kendall  Perkins 
Vice  President  (Engineering) 
McDonnell  Aircraft  Corporation 
Lambert  Street  -  St  Louis 

Municipal  Airport,  Box  515 
St.  Louis  3,  Missouri  1 

Dr.  Beno  Sternlicht 
Mechanical  Technology,  Inc, 

1  Herbert  Drive 

Latham,  New  York  3 

Mr.  J,  W.  Lower.  Chief  Engineer 
Inertial  Components 
Honeywell  Aero  Division 
2600  Ridgway  Road 

Minneapolis,  Minnesota  1 

Mr.  Carl  F,  Graesser,  Jr. 

Director  of  Research 

New  Hampshire  Ball  Bearings,  Inc. 

Peterborough,  New  Hampshire  1 

Mrs.  Alice  Ward,  Librarian 
Norden  Division  of 

United  Aircraft  Corporation 
Helen  Street 

Norwalk,  Connecticut  1 

Northrop  Corporation 
Norair  Division 
1001  East  Broadway 
Hawthorne,  California 
Attn:  Technical  Information,  1 
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Mr.  N,  L.  Simmons 
Nortronics 

A  Division  of  Northrop  Corp, 

500  East  Orangethorpe  Avenue 
Anaheim,  California  1 

Nortronics 

A  Division  of  Northrpp  Corp. 

100  Morse  Street 

Norwood,  Massachusetts 

Attn:  Mr.  E.L.  Swainson , Tech. Asst . 

Precision  Products  Dept.  1 

Pratt  &  Whitney  Aircraft 

Division  of  UAC  -  CANEL 

P.O.  Box  611 

Middletown,  Connecticut 

Attn:  Librarian  1 

Library,  Bldg.  10-2-5 
Radio  Corporation  of  America 


Camden  2,  New  Jersey  1 

Mr.  Robert  S.  Siegler 
Rocketdyne 

Nucleonics  Subdivision 
6633  Canoga  Avenue 

Canoga  Park,  California  2 

Ryan  Aeronautical  Company 
Attn:  Engineering  Library 

Lindbergh  Field 

San  Diego  12,  California  1 

Mr.  M.  A.  Vreeland 
Sanderson  &  Porter 
72  Wall  Street 

New  York  5,  New  York  1 

Jack  &.  Heintz,  A  Division  of 
The  Siegler  Corporation 
1725  Eye  Street,  Suite  505 
Washington  6,  D.C.  1 


Mr.  Paul  A.  Pitt,  Vice  President, 
Engineering  and  Research 
Solar  Aircraft  Company 
2200  Pacific  Highway 

San  Diego  12,  California  1 
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Mr.  W.  G.  Wing 

Sperry  Gyroscope  Company 

C-2 

Great  Neck,  New  York  1 

Mrs.  Florence  Turnbull 
Engineering  Librarian 
Sperry  Gyroscope  Company 
Greatneck,  New  York  1 

Sunstrand  Aviation-Denver 

2480  West  70th  Avenue 

Denver  21,  Colorado  1 

Mr.  James  W.  Salassi,  Pres. 
Turbocraft,  Inc. 

1946  S.  Myrtle  Avenue 

Monrovia,  California  1 

Universal  Match  Corporation 
Avionics  Dept.  Tech.  Library 
4407  Cook  Avenue 

St.  Louis  13,  Missouri  1 

Waukesha  Bearings  Corporation 

P.0.  Box  3^6 

Waukesha,  Wisconsin 

Attn:  Mr.  J.M.  Gruber,  Ch.Engr.  1 

Mr.  John  Boyd 

Westinghouse  Electric  Corp, 

Research  Laboratories 

East  Pittsburgh,  Pennsylvania  1 

Mr.  H,  Walter 
Director  of  Research 
Worthington  Corporation 
Harrison,  New  Jersey  1 

Dr.  W.  A.  Gross 
Ampex  Corporation 
934  Charter  Street 

Redwood  City,  California  2 

Professor  J.  Modrey 

Dept,  of  Mechanical  Engineering 

Union  College 

Schenectady  8,  New  York  2 

Stratos  Division 

Fairchild.  Stratos  Corporation 

Bayshore,  L.I.,  New  York 

Attn:  Mr.  John  Meacher  1 


A.  C.  Spark  Plug  Division 
General  Motors  Corporation 
Route  128 

Wakefield,  Massachusetts 

Attn:  Technical  Library  1 

The  Cleveland  Graphite 
Bronze  Company 
17000  St.  Clair  Avenue 
Cleveland  10,  Ohio 

Attn:  Mr.  R.  H.  Josephson  1 

Mr.  Richard  J.  Matt,  Manager 

Bearing  Development  &  Contract 

New  Departure  Division 

General  Motors  Corporation 

Bristol,  Connecticut  1 

Mr.  James  J.  Barker 
10  Walden  Avenue 

Jericho,  New  York  1 

Thomson  Ramo  Wooldridge 
TAPCO  Group 

New  Devices  Laboratories 
7209  Platt  Avenue 
Cleveland  4 ,  Ohio 

Attn:  Mr.  0.  Decker  1 

Dr.  L.  Licht 

International  Business  Machines 
Corporation 

Thomas  J.  Watson  Research  Center 
P.0,  Box  218 

Yorktown  Heights,  New  York  1 

Mr.  E.  A.  Babley 

Union  Carbide  Nuclear  Company 

P.0.  Box  P 

Oak  Ridge,  Tennessee  1 

Lockheed  Aircraft  Corporation 
Missiles  and  Space  Division 
Technical  Information  Center 
3251  Hanover  Street 

Palo  Alto,  California  1 

Office  of  Technical  Services 

Department  of  Commerce 

Washington  25,  D.  C.  1 
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